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We present the first results of an ongoing effort to determine the lattice scale on the Nf = 2 CLS 
lattice ensembles via the mass of the Q. baryon. Results from different methods are compared, 
and various sources of systematic uncertainty are discussed. 
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1. Introduction 

In order to give absolute values for dimensionful quantities (including masses and decay con- 
stants) from a lattice simulation, we need to know the lattice spacing a. 

Since no unique prescription for scale-setting exists, there is a certain ambiguity arising from 
the choice of the dimensionful quantity used to set the scale, which however is at least formally 
small since it is of higher order in the lattice spacing. Of more practical importance are the uncer- 
tainties that arise from an inaccurate knowledge of the scale caused by statistical and systematic 
errors on the relevant quantity. 

Here, we present the result of our effort to set the scale via the mass of the Q. baryon [1]. The H 
mass has a number of advantages over other quantities, in that the Q. is stable in QCD (as opposed 
to e.g. the K* meson), and that its mass is only very weakly dependent on the light quark mass (as 
opposed to e.g. the nucleon). In addition, all calculations necessary to determine the scale from 
the mass can be done in the fully relativistic theory (as opposed to e.g. a scale determination 
from the T 2S-1S splitting), and no renormalization constant is needed for ma (as opposed to e.g. 
/k)- Possible disadvantages of using the Q. baryon include the less favourable signal-to-noise ratio 
for baryons, which may render the extraction of a precise value for amQ difficult. In an A^f = 2 
simulation, an additional potential downside of the ^l baryon is the quenching of the strange quark, 
which may introduce additional (presumably small) systematic effects. 

The n baryon thus provides a way to set the scale that reduces many sources of systematic 
error in exchange for an increase in statistical error. To the extent that this increased error can be 
beaten by performing more measurements, this is a rather favourable situation. 

1.1 The CLS ensembles 

Coordinated Lattice Simulations (CLS) is a consortium designed to pool the human and com- 
puter resources of several teams in Europe interested in lattice QCD. CLS member teams are lo- 
cated at CERN, in Germany (Berlin, DESY/Zeuthen, Mainz, Wuppertal), Italy (Rome, Milan) and 
Spain (Madrid, Valencia). All CLS simulations use either M. Liischer's implementation of the DD- 
HMC algorithm [2] or the MP-HMC algorithm by M. Marinkovic and S. Schaefer [3] to efficiently 
simulate Nf = 2 Wilson QCD with non-perturbative 0(a) improvement on a variety of computer 
architectures ranging from PC clusters to the BlueGene/P at NIC/Forschungszentrum Jiilich. 

The ensembles used in the present analysis are listed in table 1 . 

2. Methods 

2.1 Measurements 

We measure the two-point correlators for mesons and baryons using Gaussian-smeaied sources 
[4] with HYP-smeared links [5] in the covariant Laplacian to suppress excited state contributions. 
The parameters of the smearing are fixed at each jS so as to correspond to a smearing radius of 
approximately 0.4 fm. Multiple point sources per configuration are used in order to improve statis- 
tics. For the vector mesons and baryons, we also average over multiple polarizations to further 
reduce statistical eiTors. 
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Figure 1: Examples of effective mass plots with plateaux and two-state fits shown for a pseudoscalar and 
vector meson and a decuplet baryon channel, corresponding to a trial K, K* and Q. on the N4 ensemble, 
respectively. 



From each measured correlator we extract the effective masses and apply the following meth- 
ods to extract estimates for the mass of the ground state: First, a naive plateau fit is performed for 
each correlator, starting at a reasonably large Euclidean time; secondly, taking the plateau value as 
the starting point, a two-state fit including the leading excited state contribution is performed with a 
larger time range. Finally, another two-state fit with the same time range is performed for the bary- 
onic channel, in which the gap between the ground and excited state is fixed to the theoretically 
expected gap of Imj^ using the measured pion mass. Examples of effective mass plots are shown 
in fig. 1. 

2.2 Fixing JCstrange 

In order to establish the physical value of the baryon, we need to fix the strange quark mass 
to its physical value. Since we intend to set the scale from the Q. mass, a scale-free renormalization 
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Table 1: CLS ensembles used in this analysis. A^^fg is the number of configurations used for the calculation 
of hadron masses in this study, not the total number generated. 
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Figure 2: An example of the interpolation to the physical value of the strange quark mass. The measured 
ratios are shown as blue diamonds, and the interpolated value as the red circle. 

condition is needed for this purpose, and we therefore will fix JCgtrange by demanding that a ratio 
of masses equals its physical value. Different choices of the ratio used correspond to different 
renormalization conditions, which will introduce an ambiguity in the scale through the resulting 
ambiguity in JCstrange- 

The three ratios which we have chosen to fix K"strange are 

• Ri = mK/niK*, which is the ratio originally used by the CERN group in [6], 

• a chirally improved variant R2 = {m\ — ^m^)/m^t, from which the leading dependence 
of Ri has been removed, and 

• /?3 = (m^ — 2m|)/m^, which is obtained by replacing the K* resonance by the stable Q. 
baryon in /?2- 

Determining all three ratios /?,• requires us to measure mj^, wip^, /My and mf^ for a range of trial 
K'strange at cach valuc of j8 and K"sea; we then determine the value of K"strange corresponding to the 
physical strange quark mass by interpolating to the physical value of An example is shown in 
fig. 2. 

The difference in the value of K^tange that results from the different mass fits is generally less 
than the statistical errors, which implies that the systematic eiTors from excited states are well 
under control. Similarly, the difference between a linear and a quadratic interpolation in fCstiange or 
in JCj^i-ange is much smaller than the statistical errors. All statistical errors were determined using 
the UWerr procedure [7] for eiTor propagation including the effects of autocoiTclations. In the 
following, we have chosen the results from the fit with a fixed gap of Irriji as our final value. 
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Figure 3: The values of fCstrange from the three renormalization conditions studied. From left to right, the 
j3 = 5.2, 5.3 and 5.5 ensembles are shown with the values obtained from Ri, R2 and Rt, plotted as blue 
diamonds, pink circles and red squares, respectively. 

On the other hand, the ambiguity in JCstrange from the choice of the renormahzation condition is 
significant, but decreases as expected towards the chiral limit, as can be seen in fig. 3. We find that 
condition 3, which employs the Q. baryon and a chirally improved estimate of the strange quark 
mass, is the most stable and therefore we choose it for our final analysis. 

2.3 Fixing the scale 

The value ro/a of the Sommer scale in lattice units has been determined on the CLS config- 
urations by Knechtli and Leder [9]. On each ensemble, we combine their results for ro/a and our 
results for am^ using the UWeiT procedure [7] in order to take into account both the correlations 
between observables and the autocoiTclations within each ensemble (for an in-depth discussion of 
the issues associated with the latter cf. [8]). 

Since our respective measurements for the jS = 5.5 ensembles were performed on non-overlapping 
subsets, we combine the results on these ensembles neglecting the correlations between ro and niQ; 
since we find that the error of rQiriQ is dominated by the error on ama, any resulting systematic 
error may be assumed to be small. 

The finite combination rotriQ can now be extrapolated to the physical point in and in 
{mj^/maY. To estimate systematic uncertainties, we used four different fits: 

• a linear fit in {mji/mo)^, assuming no lattice artifacts, 

• a quadratic polynomial in {mjt/mo)^, assuming no lattice artifacts, 

• a linear fit in {mji/mo)^ at each j3 separately, followed by a linear fit of the chirally extrapo- 
lated results in (a/ro)^, and 

• a combined continuum and chiral extrapolation by a fit linear in both {niji/mQif - and (a/ro)^. 

The first two fits, which are made under the assumption that lattice artifacts are absent from the 
combination roniQ, are motivated by the scaling behaviour seen in fig. 4. All fits agree within their 
respective errors. 
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Figure 4: Left: the joint linear and quadratic chiral extrapolations of roniQ ignoring discretization effects; 
right: the separate chiral extrapolations of roma at each value of j3. The j3 = 5.2, 5.3 and 5.5 ensembles are 
shown as blue diamonds, green circles and red squares, respectively. 



3. Results 

A combined continuum and chiral extrapolation of roniQ gives a physical value of roniQ = 
3.99(12) (9), where the errors are statistical and systematic, with the systematic error estimated 
from the spread of the results obtained using the different fitting strategies both for the mass ex- 
traction and in the extrapolation. 

Combining this with the experimental value mq = 1672.45(29) MeV [10] gives a value of 
ro = 0.47 1(14) (10) fm for the Sommer scale, which is compatible with the ETMC result [11] of 
ro = 0.465(6)(15) using the nucleon. 

Performing a separate chiral extrapolation of aniQ to the physical pion mass at each value of 
j3, we find the lattice spacings at the physical point to be 

~^ 5^ 53 5^5 

a [fm] 0.079(3)(2) 0.063(2)(2) 0.050(2)(2) 

We note that the lattice spacings obtained from the baryon mass are significantly finer than 
those previously obtained [12] for the same ensembles using the CERN method [6]. This is not 
entirely unexpected, since the CERN method uses an unphysical reference point of niji/niK = 0.85 
to make contact with the physical hadron masses, and was intended primarily for determining the 
relative scale between different j8 values rather than an absolute value of the lattice scale. 

The new values presented here are, however, compatible with the results obtained in an entirely 
separate approach from /k [13], as well as with the results from an independent analysis of the 
same data using a different eiTor estimator and performing the chiral extrapolation in [mji/mi^Y 
[14], which gives us confidence in the robustness of the results. 
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4. Conclusions 

The mass of the Q. baryon provides a way to set the scale in a way that is particularly insensitive 
to the sea quark mass, and independent of renormalization constants and effective theories. We 
find that the relatively large statistical errors of the D. baryon mass can be controlled, and that 
the systematic errors due to excited states are also well under control; hence the H provides an 
attractive means to set the scale. 

Using mci, we have determined the Sommer scale to be ro = 0.471(14)(10) fm and have set 
the lattice spacing on the Nf = 2 CLS lattices. The lattice spacings determined in this way are 
compatible with those found in an independent analysis employing /k to set the scale. We find 
that the lattices are finer than had been previously expected. What influence this will have on the 
question of finite- volume effects on studies of hadronic structure [15] remains to be seen. 

In the immediate future, we intend to increase the number of measurements on the 07 ensem- 
ble, which will reduce the statistical error on the lattice spacing at j3 = 5.5; in the longer run, we 
will also include in our analysis the G8 and N6 ensembles cun^ently being generated, which should 
help us to improve our control of the chiral extrapolation. 
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